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We produce YbF molecules with a density of 10 18 m -3 using laser ablation inside a cryogenically- 
cooled cell filled with a helium buffer gas. Using absorption imaging and absorption spectroscopy we 
study the formation, diffusion, thermalization and optical pumping of the molecules. The absorption 
images show an initial rapid expansion of molecules away from the ablation target followed by a 
much slower diffusion to the cell walls. We study how the time constant for diffusion depends on the 
helium density and temperature, and obtain values for the YbF-He diffusion cross- section at two 
different temperatures. We measure the translational and rotational temperatures of the molecules 
as a function of time since formation, obtain the characteristic time constant for the molecules to 
thermalize with the cell walls, and elucidate the process responsible for limiting this thermalization 
rate. Finally, we make a detailed study of how the absorption of the probe laser saturates as its 
intensity increases, showing that the saturation intensity is proportional to the helium density. We 
use this to estimate collision rates and the density of molecules in the cell. 



I. INTRODUCTION 

Cold molecules are useful for a wide variety of ap- 
plications in physics and chemistry. A review of these 
applications and the current status of the field is given 
in reference pQ. An important application is the use of 
cold molecules for testing fundamental laws of physics 
[2]. Cold molecules are being used to test the standard 
model of particle physics and its extensions via a mea- 
surement of the electron's electric dipole moment [3H5]. 
They can also be used to test the time invariance of the 
fine structure constant and of the proton-to-electron mass 
ratio [6HTU], test Lorentz invariance [TTj, and measure 
parity- violating interactions in nuclei [12 j and in chiral 
molecules [T3] . 

Buffer gas cooling is a very versatile way of producing 
high densities of atoms and molecules at temperatures 
down to a few hundred milliKelvin [14]. The molecules 
are introduced into a cryogenically-cooled cell where they 
thermalize with a cold buffer gas. The buffer gas is usu- 
ally helium, though neon has also been used [15] . The 
molecules can be brought into the cell in a number of 
ways, depending on their vapour pressure and reactivity. 
Methods that have been demonstrated include laser ab- 
lation of a precursor target e.g. [16], injection through 
a capillary e.g. [17], loading from a molecular beam e.g. 
[18], and direct gas flow loading e.g. [T5] . 

Once loaded, there are three main options available. 
In the simplest case the molecules are allowed to diffuse 
freely to the cell walls, where they will stick, and mea- 
surements are made on the molecules in the interval of 
time between loading and sticking. This method has been 
used to measure molecular parameters by millimetre and 
submillimetre spectroscopic techniques [T7] [I9j [20] , and 
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more recently by laser spectroscopy, e.g. [21, 22 j. Elastic 
and inelastic collision cross-sections can also be measured 
this way, e.g. [23j[24]. A second approach is to trap the 
molecules magnetically by placing the cell inside a strong 
magnetic trapping field and pumping away the helium as 
soon as the molecules are cold [25] . The coldest molecules 
are then confined in the magnetic trap. Several molecu- 
lar species have been trapped this way, [26j [27] , allowing 
elastic collision and Zeeman relaxation rates to be mea- 
sured at low temperatures, [28]. Trapping times of 20 s 
have been demonstrated for NH molecules [29 . A third 
approach is to extract the molecules from the cell to make 
a high intensity cold beam [30 which could then be used 
in a molecular beam experiment or for trapping. The 
molecules that form the beam can be guided away from 
the cell and separated from the helium using a magnetic 
guide [31], an electrostatic guide [32 or an alternating 
gradient electric guide [33] . 



It is important to have a good understanding of the 
processes that occur inside the buffer gas cell, since this 
is the starting point for all experiments whether they be 
done directly inside the cell, in a magnetic trap loaded 
from the cell, or in a beam extracted from the cell. In 
this paper we focus on the formation, diffusion and ther- 
malization of YbF molecules produced by laser ablation 
of a target inside a closed helium buffer gas cell. These 
molecules are of particular interest for the measurement 
of the electron's electric dipole moment [3 r 4 and, for that 
application, the formation of a high density of YbF is one 
of our aims. The work presented here will be applicable 
to many other buffer gas cooling experiments, particu- 
larly those where laser ablation is the loading method or 
where laser spectroscopy is used to make measurements 
inside the cell. 
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FIG. 1: (Color online) An illustration of the cell showing the 
position of the target and the gas inlet and the propagation 
directions of the ablation and probe lasers. The coordinate 
system we use is also shown and has its origin at the centre 
of the cell. 



II. EXPERIMENTAL DETAILS 

Figure [I] shows a few details of the experiment. To 
make the cell we took a solid aluminium cube of side 
40 mm and bored three orthogonal 20 mm diameter holes, 
two passing all the way through the cell and the other 
half way through. The face without a hole mounts onto 
the cold plate of a liquid helium cryostat (Infrared Lab- 
oratories, model HD3). Vacuum-tight indium seals close 
the five open faces, three with windows and two with 
blanks. One window is used for laser ablation of a target 
placed on the opposite side of the cell and the other two 
are used for laser absorption spectroscopy and imaging 
of the molecules formed inside. The cryostat is pumped 
by a turbomolecular pump to a pressure of approximately 
10 -6 mbar. A radiation shield, in thermal contact with a 
liquid nitrogen dewar, surrounds the setup and prevents 
room temperature radiation from reaching the cold sur- 
faces. Windows in this shield, coated with indium tin ox- 
ide, provide optical access at visible wavelengths, whilst 
reflecting thermal radiation. 

YbF molecules are created in the cell by laser ablation 
of a vacuum hot-pressed target consisting of 70% Yb and 
30% AIF3 by mass. The ablation pulses have a duration 
of 8 ns, an energy of approximately 50 mJ, a wavelength 
of 1064 nm, and a spot size at the target of 2 mm. 

Helium gas enters the cell through a hole beneath the 
target, via a thin- walled tube of inner diameter 2.5 mm. 
This is in two parts: a stainless steel tube for most of the 
length is connected to a copper tube near the cold plate 
of the cryostat. To ensure that the gas is cold when it 
enters the cell, the copper tubing winds around, and is 
silver-soldered to, a copper bobbin connected to the cold 
plate. At the room temperature end, outside the cryo- 
stat, the feed tube is connected to a chamber where the 
helium pressure is measured and controlled. The helium 
enters this chamber through a leak valve and is pumped 




Inlet pressure (mbar) 

FIG. 2: (Color online) Pressure measured at the cell versus 
pressure measured at the inlet tube, all at room temperature. 
Red squares: during ablation. Black open circles: pressure 
settled after ablation. Blue filled circles: after pumping on 
the cell through the gas inlet tube for about 12 hours. The 
errors are about 1% of the reading. The line has zero offset 
and a gradient of 1. Note the logarithmic scale. 



away to a rotary pump through a second valve. These 
two valves provide fine control over the helium pressure, 
which is measured by a Pirani gauge. When the cell is 
cold, there is a thermomolecular pressure drop between 
the two ends of the tube, and we account for this using 
the empirical formula given in [34 . Even when the cell is 
at room temperature, it is not immediately clear whether 
the pressure measured at the inlet end of the feed tube is 
the same as the pressure in the cell. To investigate this, 
we temporarily connected a second Pirani gauge directly 
to the cell, via a port in the cryostat usually used as a 
window. Both pressure gauges were calibrated to a third 
gauge which, for helium, is accurate to 0.5%. 

Figure [2] shows the pressure in the cell versus the pres- 
sure at the inlet. One set of measurements was taken 
during ablation of the target at 10 Hz, another set a few 
minutes after stopping the ablation, and a final set after 
pumping the cell through the gas feed tube for about 12 
hours. The figure shows that the cell pressure is equal to 
the inlet pressure plus a constant offset. The offset is due 
to outgassing of the target which raises the total pressure 
in the cell. This outgassing is greatest during ablation of 
the target, when the offset is 0.17 ± 0.05 mbar. When 
the ablation is stopped the offset falls with two sepa- 
rate timescales. First it falls quickly, reaching a steady 
value of 0.06 ± 0.02 mbar in a few minutes, and then it 
falls much more slowly, reaching 0.015 ± 0.003 mbar af- 
ter about 12 hours. We could only make these direct 
pressure measurements at room temperature, so we do 
not know how the offset changes with temperature. The 
pressure offset has implications for some of our measure- 
ments, and we will highlight these where relevant. 

The temperature is measured using a silicon diode tem- 
perature sensor attached to the cold plate of the cryostat. 
At the lowest temperatures however, temperature mea- 
surements of the fully thermalized molecules (Sec.jv]) in- 
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dicate that the cell is about 10 K hotter than the cold 
plate, presumably due to insufficient thermal conduction 
between the two. We have also used the Doppler broad- 
ening of the atomic Yb spectrum to determine the tem- 
perature, and this agrees with the results obtained from 
the molecules. Using a new experimental setup, where 
the thermal conduction is improved, we find that the 
molecules do indeed reach the same temperature as the 
4K cell. All the measurements reported here were done 
with the earlier apparatus and we use the temperature 
of the molecules as the most reliable measure of the cell 
temperature. 

All the measurements presented in this paper are based 
on absorption of a cw probe laser beam. Light from a 
dye laser is used to address individual rotational lines 
within the (0-0) band of the X 2 S + -A 2 n 1 / 2 transition 
of YbF, at a wavelength of 552 nm. The laser frequency 
is measured with an accuracy of approximately 600 MHz 
using a wavemeter (HighFinesse WS6). Absorption im- 
ages reveal the spatial distribution of the molecules as 
a function of time. For these measurements, the laser 
beam passes through a 100 MHz acousto-optic modulator 
(AOM) acting as a fast optical shutter, and is then colli- 
mated and expanded to a diameter of 20 mm. This beam 
passes through the buffer gas cell and onto a ccd camera 
(Marlin F033B). The shortest exposure time offered by 
this camera is 34 /is. To take images with higher time 
resolution, we use the AOM to turn on the absorption 
imaging beam for a period of 10 /is. For time-resolved 
measurements of the YbF density, translational temper- 
ature and rotational temperature at a specific point in 
the cell, we use a combination of Doppler-limited and 
Doppler-free absorption spectroscopy, using 1 mm diam- 
eter pump, probe and reference beams, as described in 
detail in [22]. 

Data are taken at a repetition rate of 10Hz, synchro- 
nised to the 50Hz line frequency to suppress the effects 
of line noise in the experiment. In the imaging experi- 
ments the laser frequency is fixed and each shot yields 
one image from the ccd camera. In the spectroscopy ex- 
periments, the outputs of the photodiodes are recorded 
for 8 ms with a sample rate of 250 kHz, and the laser fre- 
quency is stepped between one shot and the next. As de- 
scribed in [22] , the absorption of a fixed frequency probe 
beam from a second dye laser is used for normalization 
purposes, factoring out the slow drift in ablation yield 
that occurs as the target degrades. 



III. ABSORPTION IMAGING 

Figure [3] shows a sequence of absorption images for 
three different buffer gas densities at a cell temperature 
of 20 K. To obtain each image, we take the difference 
between two pictures, one where the ablation laser fires 
and one where it does not. The field of view is 20 mm 
wide and the colour represents the fractional absorption 
of the probe laser. The target is on the left side of the 



image, just outside the field of view. The probe laser 
is tuned near the P-branch bandhead where there is a 
high density of overlapping spectral lines from rotational 
states in the range N = 3 — 6. The camera is exposed to 
the probe for 10 /is with a variable delay relative to the Q- 
switch of the ablation laser. The images therefore show 
how the density of YbF molecules in low-lying rotational 
states evolves with time, and how this depends on the 
buffer gas density. 

In Fig.[3^a), the helium density is relatively low, nn e = 
0.9 x 10 25 m- 3 . In this case, the molecules fill the entire 
field of view within 30 /is. Their distribution is roughly 
uniform in all the images. The overall absorption in- 
creases for the first 180 /is, and then slowly decreases 
after that. We propose that the molecules are formed 
close to the target, then expand ballistically, and because 
their mean free path is relatively long at this density, they 
rapidly fill the entire cell. They then diffuse slowly to the 
walls, where they stick. The increase in absorption over 
the first 180/is suggests either that the formation of new 
molecules continues for this length of time, or that the 
rotational temperature starts out high and cools on this 
timescale, thereby increasing the population of the low- 
lying rotational states being probed in the experiment. 
The measurements presented in Sec.[V] confirm that the 
rotational temperature does indeed fall considerably over 
this initial period. In the first frame of the sequence of 
images hardly any molecules are visible because they are 
too hot during this initial 10 /is period to be observed. In 
Fig.[3|b) the helium density is increased to 17 x 10 22 m -3 , 
and the evolution has changed considerably. The initial 
ballistic expansion is arrested before the molecules have 
filled the cell, so that in the first frame we see the pro- 
jection of a hemispherical distribution with a radius of 
about 8 mm. This distribution then diffuses slowly into 
the rest of the cell. The absorption peaks after about 
60 /is, indicating that the rotational cooling is faster than 
before. We might expect the decay of the molecule den- 
sity to be slower than before because the diffusion will be 
slower at higher helium density, but in fact the density 
appears to decay more rapidly. This is even more evi- 
dent in Fig.[3tc) where the helium density is increased to 
94 x 10 22 m". Here, the ballistic expansion is arrested 
even closer to the target, the peak absorption occurs at 
even earlier times, and the propagation into the rest of 
the cell is even slower. After 540 /is the molecules have 
still not filled the entire field of view, and yet most of the 
molecules have vanished. The rapid decay of the signal 
is partly due to the proximity of the molecules to the 
walls of the cell. The target is close to the wall, and the 
molecules are stopped close to the target, so the distance 
that they have to diffuse before they reach a wall is small. 
This tends to lower their survival time in the cell, even 
though the mean free path is small. Other loss processes 
may also be contributing to the decay of the molecular 
density. The YbF molecules and other ablation products 
are confined at high density and the YbF radicals may 
react with these other atoms and molecules to form more 
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(b) 17 x 10 22 m- 3 



(c) 94 x 10 22 m 3 





360 jas 



450 jas 



540 jas 



FIG. 3: (Color online) Absorption images of YbF at various times since ablation and for three different helium densities at a 
cell temperature of 20 K. 



stable products. We investigate the decay times of the 
molecules more quantitatively in the next section. 



IV. DIFFUSION 

We have studied the diffusion of the YbF molecules 
through the helium buffer gas. In order to obtain an 
estimate of the YbF-He diffusion cross-section, without 
needing an accurate measure of the helium density in 
the cell, we produce and probe lithium atoms in the cell 
simultaneously with the YbF molecules. The Li-He dif- 
fusion cross-section can be calculated accurately, and so 
the Li diffusion times provide a good reference against 
which the YbF diffusion times can be compared. 



where /i = mM / (m + M) is the reduced mass related to 
the masses of the helium, m, and the diffusing molecules, 
M. Although this expression is only the first term of 
a series approximation for the diffusion coefficient, it is 
expected to be accurate to better than 2% for the two 
cases, Li-He and YbF-He, that we consider here [37] . 



The density of molecules in their journey through time 
and space is described by the time dependent diffusion 
equation ^ = \7 2 (Dn). The solution can be written in 
the form n(r, t) = Ckfk(r)e~ t / Tk , where is the am- 
plitude of the diffusion mode /fc(r), which for a homoge- 
neous diffusion coefficient is a solution of the differential 
equation 



A. Theory 

Following ablation, there is an initial ballistic expan- 
sion of molecules away from the target. Over a certain 
distance from the target, depending inversely on the he- 
lium pressures, this ballistic expansion is arrested and 
the molecules then diffuse towards the walls. The flux 
of diffusing particles, J, is proportional to the gradient 
of the density n, J = DVn. For diffusion of one gas 
into another at temperature T, the relationship between 
the diffusion coefficient, D, and the thermally averaged 
diffusion cross-section, <td, is given, with sufficient accu- 
racy, by the first Chapman-Enskog approximation [35]- 
[37] . Since the number density of helium, nn e , is many 
orders of magnitude larger than that of our diffusing gas, 
the result is 

16cr D n H e V /i 



y^+V 2 A = (2) 



subject to the boundary conditions imposed by the walls. 
Since the molecules stick to the walls with high proba- 
bility, we take n = at the walls. When there is no 
particular symmetry in the problem, three indices are 
needed to label the diffusion modes, represented here by 
the single label k. In all cases, the product of the diffusion 
coefficient and the time constant depends only on the 
cell geometry and on the indices labelling the diffusion 
mode. The higher order diffusion modes have smaller 
time constants. 

If the interaction potential between the two colliding 
species is known, the thermally-averaged diffusion cross- 
section (Td can be calculated. For all temperatures, T, 
of interest here, the calculation can be done classically 
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using the following set of equations [37] : 
1 f°° 

a B = - x 2 e~ x a B (E)dx, 
z Jo 

poo 

a D (E)=27r / [1 - cos x(E,b)]bdb, 
Jo 



(a) |as 



X (E,b)=ir-2b 



poo 



r dr 



y/l-V(r)/E -b 2 /r 2 ' 



(3) 
(4) 
(5) 



Here, <jb(E) is the diffusion cross-section at centre-of- 
mass collision energy E, x = E/(kT), b) is the 

deflection angle for a collision with energy E and im- 
pact parameter 6, V(r) is the interaction potential as 
a function of the particle separation r, and r c is the 
distance of closest approach in the collision, given by 
1 - V(r c )/E - b 2 /r 2 c = 0. 

The interaction potential for Li-He is well established, 
there being a long history of theoretical and experimen- 
tal work on this system (see [38] and references therein). 
We have calculated od, Li-He for a range of temperatures 
using the potential given in reference [38] . The results at 
293, 80 and 20 K are given in Table [I] on page [9] There 
are two sources of error in this calculation. The first is 
due to our neglect of quantum effects. Using the tabu- 
lation in [37 , we estimate that quantum effects increase 
the cross-section by about 2% at 20 K, and are negligible 
at the higher temperatures. The second source of error 
is due to the uncertainty in the interaction potential. To 
estimate this we repeated the calculations using the alter- 
native potential given in [39] . The cross-sections obtained 
from this latter potential are smaller by 0.3% at 293 K, 
5% at 80 K and 6% at 20 K. Taking these differences as 
indicative of the likely accuracy, we assign a fractional 
uncertainty of 5% to the cross-sections calculated at all 
temperatures. The interaction potential for YbF-He has 
also been calculated [40] and we use this to obtain theo- 
retical values for <jd, YbF-He at these same temperatures. 
In this case, the potential is not only a function of the 
YbF-He separation but also of the angle between the YbF 
internuclear axis and the incident velocity vector of the 
collision. We calculate an approximate diffusion cross- 
section by calculating the cross-section as a function of 
this angle and then averaging over all angles. The results 
are given in Table [T] 



B. Model 

Since our cell geometry is not simple, we model the dif- 
fusion using finite element software 1 . This model uses the 
same cell geometry as the experiment and only accounts 
for diffusion. It does not include the initial momentum 
away from the target that the molecules have following 




FIG. 4: (Color online) Simulated absorption images for a 
diffusion coefficient of 0.005 m 2 s _1 . This corresponds, for ex- 
ample, to a cell temperature of 20 K, a diffusion cross-section 



of 80 x 10 nr and a helium density of 3 x 10 m 



i Uni- 



1 Comsol Multiphysics 3.2 



form initial distribution throughout the cell, (b) Distribution 
localized near the target. 



ablation. For low helium densities, this initial momen- 
tum helps to distribute the molecules throughout the cell 
faster than diffusion would be able to do. To understand 
the effect of this, we study two extreme cases for the ini- 
tial condition in the model. In one case, the molecules 
start out with a uniform distribution over the whole cell. 
This is intended to model the situation at low helium 
density where the initial ballistic expansion completely 
fills the cell with molecules. In the other extreme case 
the initial density is taken to be e~ r l w where r is the 
distance away from the target and w = 1 mm. This is 
intended to model the situation at high helium density 
where the ballistic expansion is rapidly arrested with the 
molecules localized close to the target. 

Figure [4] shows the predictions of these simulations for 
the two extreme initial conditions. The behaviour seen 
in these simulated images is qualitatively similar to that 
of the experimental absorption images shown in Fig. [3] 
In Fig.[4ja) the initial density distribution is a constant. 
This evolves to a distribution with a maximum at the 
middle, and that distribution then decays away slowly. 
In FigQb) the initial density is localized close to the 
target, and we see that as time progresses the distribution 
expands slowly away from the target but decays away on 
a faster timescale than in case (a). This is because the 
molecules are localized close to the target and so can 
diffuse back onto the target and the nearby wall. The 
diffusion equation tells us that high spatial gradients in 
density decay away rapidly. The highly localized initial 
distribution is represented by high-order diffusion modes 
localized close to the target, and these decay away quickly 
with a large fraction of the molecules quickly reaching 
the target and nearby wall, leaving only low-order modes 
whose amplitudes are so small that they are barely visible 
in the images. 

As well as taking absorption images, we measure the 
absorption locally using a small probe laser. To simulate 
these measurements, we integrate the simulated density 
over the volume of the probe laser. Doing this, we find 
that the time evolution of the density is sensitive to the 
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position of the probe, and to the initial condition, as 
shown in Fig. [5] In this figure, we plot the simulated 
time evolution of the density probed at three different 
distances from the target, as well as the evolution of the 
density averaged over the whole cell. With reference to 
the coordinate system shown in Fig. [I] the probe laser 
is in the x?/-plane, displaced along y. In each case, the 
density is normalized to the average density in the cell at 
t = 0, and is plotted on a logarithmic scale. Our choice 
of time axis in this plot requires some explanation. Con- 
sider again the diffusion equation and introduce scaled 
coordinates r' related to the original coordinates r by 
y' = r/L where L is a characteristic length scale for the 
cell. Its precise value is unimportant and we will take 
L = 1 cm. Taking derivatives with respect to the scaled 
spatial coordinates the diffusion equation is = V 2 ,n 
where t' = Dt/L 2 is a dimensionless time. This is the 
time axis we use in Fig. [5] because this makes it obvi- 
ous how the results scale with the overall size of the cell 
and with the diffusion coefficient (and hence the diffusion 
cross-section, helium density and temperature). 

In the case where the initial density distribution is 
uniform [Fig.[5ja)], the time evolution does not depend 
strongly on position. There are some differences in the 
time evolution for different beam position at very early 
times, but after that the evolution is well described by 
a single exponential decay with a time constant which 
is the same at all positions. For our cell, this time con- 
stant is found to be r = 0.22 x 10 -4 /D, where D is the 
diffusion coefficient in m 2 s _1 and r is in s. For exam- 
ple, for a helium density of 10 22 m -3 , a YbF-He diffusion 
cross-section of 80 x 10 -20 m 2 and a temperature of 20 K, 
the diffusion coefficient is D = 0.015 m 2 s _1 , and the ex- 
pected diffusion lifetime is r = 1.5 ms. For diffusion in 
a cube of side a, the time constant associated with the 
lowest order diffusion mode is r = a 2 /(37r 2 D). It follows 
that the diffusion time constant for our cell is the same 
as that of a cubic cell with a = 2.6 cm. This is reassur- 
ingly similar to the size scale of our cell. The position 
sensitivity is far stronger when the molecules are initially 
concentrated near the target, as in Fig.^b). While the 
initial density is largest close to the source of molecules, 
there is a rapid decay of the density at this position. As 
time goes on, the rate of decay slows down, and even- 
tually the time evolution becomes a single exponential 
with the same time constant as in the constant density 
case. However, that does not occur until the density has 
fallen by two orders of magnitude. Further away from 
the source, the density first increases with time as the 
molecules arrive at the probe, and then it decreases again 
with the same time constant as in the constant density 
case. The peak density decreases as we move further 
from the target, and the time taken to reach this peak 
increases with increasing distance. We conclude that at 
long times the time evolution is insensitive to either the 
initial condition or the position of the probe, but that at 
short times it is sensitive to probe position particularly 
when the initial density is strongly localised. 
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FIG. 5: (Color online) Simulated evolution of the molecule 
density at different positions in the cell. The density is nor- 
malised to the initial density averaged over the whole cell. 
The time axis is expressed in the dimensionless units Dt/L 2 , 
where D is the diffusion coefficient and L — 1 cm is a char- 
acteristic size scale for the cell, (a) Uniform initial density 
throughout the cell, (b) Initial density localised near the tar- 



These observations can be explained by expressing 
the density distribution as a sum over a set of diffu- 
sion modes, the higher order modes representing more 
rapid spatial variations in density and decaying away 
more quickly. When the cell is uniformly filled, the den- 
sity rapidly settles into the lowest order diffusion mode, 
which gives us the longest diffusion time. To represent 
an initial distribution that is strongly localized we need 
to put large amplitude into high-order modes, and these 
then decay away rapidly. Eventually, only the low-order 
modes, which persist for much longer, will remain, but 
the amplitude of the low-order mode may be very small. 
In an experiment, it may be too small to observe above 
the noise. 



C. Experimental results 

We compare the YbF diffusion times to the diffusion 
time of lithium atoms introduced into the cell at the same 
time. If we assume that the initial distribution is the 
same for the YbF and the Li, the product of the diffusion 
coefficient and the time constant is the same for both 
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species, since this product is only a function of the cell 
geometry. It follows that 



V n 

where subscript 1 refers to Li- He and subscript 2 to YbF- 
He. To introduce Li into the cell, a piece of lithium 
wire was mounted into the target holder along with the 
AlFa/Yb target. The two materials were ablated simul- 
taneously using a single Nd:YAG laser spot. The lithium 
density was monitored by absorption of light from a diode 
laser tuned into resonance with the 7 Li Dl line at 671 nm. 
This probe laser was overlapped with the YbF probe 
beam and the absorption signals for the two species were 
measured simultaneously on two separate photodiodes. 
In this way, both species were measured at the same place 
and at the same time, and both were produced in the 
same place on the target. The probe lasers were in the 
#2/-plane (see Fig.[T]), about 7 mm from the target. It was 
possible to obtain extremely high optical depths for Li, 
but through judicious positioning of the ablation spot on 
the target we were able to obtain optical depths of the 
same order of magnitude for the two species. Because 
the Li and YbF densities are many orders of magnitude 
lower than the He density collisions between them are 
not important. 

Figure [6] shows the time evolution of the YbF and Li 
absorption signals for a helium density of 2.2 x 10 22 m -3 . 
At this density we expect the molecules to fill the en- 
tire cell through ballistic expansion. We plot the absorp- 
tion coefficient, a, which is directly proportional to the 
molecule density, the interaction length and the absorp- 
tion cross-section, and is related to the fractional absorp- 
tion A via a = — In (1 — A) (see appendix |A|). Transients 
due to optical pumping and collisional redistribution of 
the population are too fast to be observable in the ex- 
periment. The YbF density is probed near the A-X(O-O) 
bandhead where there is a dense cluster of overlapping 
P-branch lines. For Li, there is an initial fast rise in 
the density during the first 20 [is as the atoms fill the 
cell. After that, the signal is described well by a double 
exponential decay, one exponential having a short time 
constant that dominates during the first 150 /is, and the 
second having a longer time constant that dominates at 
later times. This same behaviour is observed over a wide 
range of buffer gas densities. We find that as the helium 
density increases the short time constant decreases from 
about 40 /is when the density is 0.25 x 10 22 m -3 to about 
20 /is when the density is 2.5 x 10 22 m -3 . We suppose 
that this fast decay is the decay of higher-order diffu- 
sion modes, the result being sensitive to the initial con- 
dition and the position of the probe laser as discussed 
above. The long time constant on the other hand in- 
creases with buffer gas density, and we take this to be 
the characteristic diffusion time for Li in the cell. The 
YbF absorption curves are influenced by both diffusion 
and by rotational cooling which, at early times, increases 
the relative population in the low-lying rotational states 
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FIG. 6: (Color online) Absorption coefficient on a logarithmic 
scale plotted as a function of time for Li (red) and YbF (blue) 
at a helium density of 2.2 x 10 22 m -3 and a cell temperature 
of 293 K. 



probed in the experiment, particularly for low cell tem- 
peratures. In Sec. [V] we measure how the rotational tem- 
perature changes with time. A model that includes both 
diffusion and the rotational cooling that we have mea- 
sured predicts that, for a cell temperature of 20 K, the 
change in absorption is strongly influenced by rotational 
cooling at early times, but is dominated by diffusion at 
later times. At higher temperatures the influence of rota- 
tional cooling is far weaker. To obtain the YbF diffusion 
time constants we fit single exponential decays to the late 
part of the absorption curves (t > 0.32 ms for the 20 K 
data and t > 0.14 ms for 293 K). Our model shows that 
we will obtain accurate diffusion times using this proce- 
dure. We have also verified that, at 20 K, fitting only to 
data beyond 0.6 ms makes a negligible difference to our 
results. 

Figure [7] shows the characteristic diffusion time con- 
stants for Li and YbF as a function of helium density, 
for cell temperatures of 293 K, 80 K and 20 K. At 80 K 
no useful YbF data was obtained, so only Li data are 
shown. The error bars are determined from the scatter 
of repeated measurements taken under the same nominal 
conditions. They therefore appear as errors in the time 
constant though we believe the major source of the scat- 
ter to be due to the difficulty of setting the buffer gas 
density, which was very slow to reach equilibrium partic- 
ularly at low temperatures. For low buffer gas densities, 
all the measured time constants increase linearly with 
density and both species diffuse more slowly at lower tem- 
peratures as expected in a simple model of diffusion. At 
both 293 K and 20 K the time constant for YbF increases 
more rapidly with increasing helium density than for Li. 
At higher densities however, the diffusion times stop in- 
creasing and become roughly independent of the helium 
density. This behaviour is displayed by both YbF and Li 
at 293 K, by Li at 80 K and YbF at 20 K. The data for Li 
at 20 K appear to show a more complicated behaviour 
at high density. For Li, the lifetimes level off at ap- 
proximately 400 /is at 293 K, 600 /is at 80 K, and roughly 
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FIG. 7: (Color online) Diffusion time constants as a function 
of buffer gas density for Li (open circles) and YbF (filled cir- 
cles) for buffer gas temperatures of (a) 293 K, (b) 80 K (Li 
only) and (c) 20 K. The lines (solid for Li, dashed for YbF) 
are linear fits to the data at low buffer gas density Note the 
different scales on both axes for the three graphs. 



two fits are consistent within errors, and we take the 
mean value as our best estimate. The gradients of the 
linear fits, together with Eq. (|6| and the calculated Li- 
He diffusion cross-section, provide a measurement of the 
YbF- He diffusion cross-section, £Td, YbF-He- The results 
are given in the fourth column of Table [IJ The uncer- 
tainty is derived by adding in quadrature the statistical 
errors in the gradients of the linear fits, the error in the 
Li-He cross-section estimated above, and in the case of 
the 20 K result, half the difference between the gradients 
of the 4-point and 5-point fits. Neither the helium den- 
sity or the cell temperature enter in the determination 
of the cross-section. Errors arising due to a difference in 
the initial spatial distribution of the Li and YbF or small 
displacements of the two probes are deemed to be negligi- 
ble because at these helium densities the initial ballistic 
expansion fills the cell, as in Fig.[3ja), and our simula- 
tion results indicate that the diffusion time is insensitive 
to the exact initial distribution or probe position in this 
case. Errors arising due to a difference in temperature be- 
tween the two species are also negligible because the time 
for thermalization with the helium is much smaller than 
the characteristic diffusion time (see Sec.|V|). Our exper- 
imental values for the YbF-He diffusion cross-section are 
consistent with our theoretical values. In particular, the 
values are in agreement at 293 K where the experimental 
uncertainty is smallest. 

Extrapolating the data to zero helium density as shown 
by the lines in Fig. [7] we observe an offset in all of the 
lifetimes. This suggests that the ablation process itself 
puts a significant density of material in the cell which 
the YbF and Li then have to diffuse through. This is 
consistent with the offset observed in Fig. [2] The offset 
is smaller for the Li data at 80 K. We do not know the 
reason for this. We speculate that, when the cell is warm, 
vacuum contaminants such as water and hydrocarbons 
reach the cell and are absorbed into the target, and when 
the cell is very cold the target absorbs helium just like a 
cryo-sorb. At 80 K neither occur because this is not cold 
enough to absorb helium but is too cold for water and 
hydrocarbons to reach the cell. 



D. Discussion 



500 /is at 20 K. For YbF, the corresponding values are 
approximately 800 /is at 293 K and 1000 /xs at 20 K. The 
density at which the levelling off occurs is about the same 
for the two species, and shifts to lower values as the tem- 
perature is lowered. 

At low density, the diffusion times show the expected 
linear increase with helium density, so we fit straight lines 
to the low density parts of the five datasets. The fit re- 
sults are shown by the lines in Fig.[7| For the 293 K and 
80 K data, the first 10 data points were used for the lin- 
ear fit. For the 20 K data it is unclear whether to fit 
to the first 4 or first 5 data points (the latter is shown 
in Fig.[7^c)), and so we do both. The gradients of these 



We have used the Li data to calibrate the YbF data 
and so our determination of the YbF-He diffusion cross- 
section does not rely on the diffusion model or on know- 
ing the helium density. It is nevertheless interesting to 
compare the experimental results at low density with the 
predictions of the diffusion model for Li. We assume that 
at low helium densities, the cell is immediately filled with 
a uniform density of particles. Fitting a sum of two ex- 
ponentials to the results of the model with this initial 
condition gives the longer of the two time constants as 
r = 0.22 x 10 _4 /D, where D is the diffusion coefficient 
in m 2 s _1 and r is in s. As noted above, the long time 
constant is insensitive to the initial condition, so the er- 
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TABLE I: Diffusion cross-sections and lifetimes for Li and YbF at three temperatures. The second and third columns give 
theoretical values for the Li- He and YbF-He diffusion cross-sections, calculated using the potentials of 38 and [40] respectively. 
The fourth column gives experimental values for the YbF-He diffusion cross-section obtained from the data shown in Fig. [7] 
calibrated by the Li-He cross- sections. The fifth and sixth columns give theoretical and experimental values for the Li diffusion 
lifetime divided by the helium density. The theoretical values are determined from the theoretical cross-sections and a numerical 
model of diffusion in the cell. The experimental values are the slopes of the straight line fits shown in Fig. [7] 



ror arising from the uncertainty in the initial distribution 
is small. Using the definition of D given in Eq. and 
the calculation of the Li-He diffusion cross-sections, we 
obtain a theoretical prediction for the gradient Tu/nn e . 
The theoretical results obtained at 293, 80 and 20 K are 
given in Table [I] where they are compared with the exper- 
imental gradients obtained from the fits shown in Fig.[7| 
The error bars given here are the fitting errors, and do 
not include a possible systematic error in the determina- 
tion of the helium density. The experimental and theo- 
retical results are in good agreement at 80 K and 20 K, 
but disagree by a factor of approximately 2 at 293 K. We 
do not have an explanation for this discrepancy. The 
temperature only enters the model through the diffusion 
coefficient, so if the discrepancy is due to an error in 
the model it must be that the calculated diffusion cross- 
section is accurate at 80 K and 20 K, but not at 293 K. 
This seems unlikely. The alternative is that there is a sys- 
tematic error in the experiment which is large at 293 K, 
but small at the lower temperatures. The most obvious 
source of systematic error is in the determination of the 
helium density, but here we would expect the error to be 
smallest at 293 K since the measurements were carefully 
calibrated at this temperature. We note that it is quite 
common to determine diffusion cross-sections from diffu- 
sion times in a buffer gas cell. To do so without using a 
reference atom requires a careful calibration of the buffer 
gas density and, unless the shape of the cell is particu- 
larly simple, a numerical model of diffusion in the actual 
cell geometry being used. 

Next, we turn to the question of why the diffusion 
times flatten off at higher helium density. A similar be- 
haviour, suggestive of an additional loss mechanism, has 
been noted in several other papers e.g. [16 , 23 , 241(41]. In 
reference [41 the decay times were measured as a func- 
tion of helium density for lithium, rubidium, silver and 
gold. In all cases, the decay times first increased linearly 
with helium density, but then turned over and gradually 
decreased with increasing density. The authors suggest 
three possible loss channels that might be responsible for 
the shortening of the lifetime at high density: the for- 
mation of dimers, capture of atoms by clusters formed 
during the ablation process, or atom loss on impurities 
present in the buffer gas or ablated off the target. One of 



these mechanisms might be responsible for the behaviour 
we observe. Here, we suggest a fourth mechanism, based 
only on diffusion and on the way the initial molecule dis- 
tribution depends on helium density. The absorption im- 
ages, Fig.[3j show that there is a rapid initial ballistic ex- 
pansion of molecules into the cell which sets the initial 
condition for the subsequent diffusion to the walls. When 
the density is low, the molecules rapidly fill the entire cell 
and so the initial distribution is set by the cell geometry 
and is independent of density. In this regime the diffu- 
sion time constant is proportional to the helium density. 
At higher densities however, an initial cloud of molecules 
is produced whose size depends on the helium density 
rather than on the cell size. The higher the density, the 
smaller the size of the initial cloud. This localization 
of the cloud results in a more rapid initial decay, as in- 
dicated in Fig.[5j particularly if the probe is moved to a 
place where the signal is initially high, which it is natural 
to do experimentally. We see that at high density there 
are two competing density-dependent effects: increasing 
the density decreases the diffusion coefficient which slows 
down the diffusion, but increasing the density also re- 
duces the initial size of the distribution which in turn 
drives faster diffusion, at least at the beginning. The 
turn-over from the proportional to the weakly-dependent 
regime occurs once the initial ballistic expansion fails to 
fill the cell, and this will depend on the cell's size and 
on the energy of the ablation pulse. Ideally, it would 
be possible to observe the multi-exponential decay that 
will occur at high density, and to find the longest decay 
time which will be independent of the initial condition. 
In practice the amplitude of the lowest order mode may 
be too small to observe, in which case the initial condi- 
tion will inevitably affect the measurement of the time 
constant in the way we have described. 



V. THERM ALIZATION 

The translational and rotational temperatures can be 
determined by recording Doppler-broadened spectra of 
the YbF molecules. The translational temperature is 
found from the Gaussian widths of the lines, and the 
rotational temperature is determined from the relative 
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strengths of transitions originating from different rota- 
tional states. Since we record the time-evolution of the 
absorption at each laser frequency, we can determine the 
translational and rotational temperatures as a function 
of time. 



A. Translational temperature 

The measurement of the width and the area of the 
Doppler-broadened lines is complicated by hyperfine 
splittings that are comparable to the Doppler widths of 
the features themselves, and by the existence of several 
YbF isotopic species. The use of Doppler-free saturated 
absorption spectroscopy enables us to identify the cen- 
tre frequencies of the underlying components, and hence 
to fit a combination of multiple Voigt profiles to the 
Doppler-broadened spectra [22] . Figure [8] shows typical 
Doppler-broadened spectra in an experiment where the 
cell temperature was approximately 20 K and the helium 
density was approximately 0.9 x 10 22 m -3 . Two spectra 
are shown, one averaged over a 200 /is wide time-window 
centred 200 /is after the ablation pulse, and the other 
averaged over a window of the same width but centred 
2000 /is after the the ablation pulse. In the early time 
window the presence of multiple hyperfine components is 
hidden by the large Doppler width, but this is revealed 
in the later time window when the translational tem- 
perature is lower. Underneath the Doppler-broadened 
profiles is shown the saturated absorption spectrum that 
was recorded simultaneously. There are four hyperfine 
components, though two of them remain unresolved in 
the Doppler-free spectrum so we only see three resolved 
spectral lines. Knowing the spacing of the hyperfine 
components from the Doppler-free spectrum, we fit to 
the Doppler-broadened spectra using a sum of Voigt pro- 
files, where the only free parameters are a single Gaus- 
sian width and the amplitudes of the various components. 
The width of the Lorentzian in the Voigt profile is fixed at 
30 MHz which is the typical width obtained from fitting 
Lorentzians to a large number of lines in the Doppler-free 
spectrum. The fitted widths give translational temper- 
atures of 76±4K and 23±1K for the two sets of data 
shown, clearly indicating the translational cooling of the 
YbF towards the temperature of the cell walls. 



B. Rotational temperature 

The rotational temperature of the molecules is deter- 
mined from the relative strengths of absorption lines orig- 
inating from different rotational states N. We obtain the 
relative line strengths from the areas of the Voigt profiles 
fitted to the Doppler-broadened spectra. Figure [9] shows 
the relative line strengths of the Q-branch transitions 
with N between 4 and 13, determined at two different 
times, 400 /as and 2000 /as after the ablation pulse. The 
width of the time window is 200 /as in both cases. The fits 
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FIG. 8: (Color online) Upper traces: Doppler-broadened ab- 
sorption spectra of the Q(7) transition with a cell temper- 
ature of approximately 20 K. Various hyperfine components 
and YbF isotopologues contribute to the absorption. The 
dark red markers correspond to a spectrum recorded 200 /is af- 
ter the ablation pulse that produces the molecules. The light 
blue markers show data recorded 2000 /is after ablation. The 
solid lines are multi-component Voigt fits to the data. The 
dashed lines isolate the components in the fit corresponding to 
the 174 YbF isotopologue. Lower trace: Saturated absorption 
spectrum recorded simultaneously allowing the line centres of 
the hyperfine components to be determined and fixed in the 
fits to the Doppler-broadened spectra. 



shown in the figure are Boltzmann distributions, modi- 
fied by the slight dependence of the transition matrix 
elements on the rotational state. These fits determine 
the rotational temperature to be 37±3 K at 400 /as and 
20±1 K at 2000 /as indicating that the rotational temper- 
ature of the YbF also cools over this timescale. 



C. Thermalization models 



Figure 10 shows the evolution of the translational and 
rotational temperatures plotted at intervals of 100 and 
200 /is, respectively. The error bars indicate the uncer- 
tainties in the parameters obtained from fits such as those 
shown in figures [8] and [9j The figure shows that the trans- 
lational and rotational degrees of freedom are in close 
equilibrium throughout the cooling process. The data 
also indicate that the molecules initially cool rapidly dur- 
ing the first 200 /is, and then continue to cool but on a 
much longer time scale. 

We fit this data to a model in which the temperature, 
T, of the molecules evolves with time according to 



T = T c J r Tue 



-t/T! 



T l2 e 



-t/r 2 



(7) 



This model fits very well to our data, as indicated by the 
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FIG. 9: (Color online) Relative strengths of Q-branch transi- 
tions originating from different rotational states N, at 400 fis 
after ablation (filled points) and 2000 /is after ablation (open 
points). The line strengths are determined from the areas 
of the Doppler-broadened spectral lines. The cell tempera- 
ture was approximately 20 K and the helium density approx- 
imately 0.9 x 10 22 m -3 . The solid and dashed lines are fits to 
the data, with the rotational temperature as a free parameter 
in the fit. 



best fit lines shown in Fig.[l0] Indeed, the quality of the 
fits to the data tend to be better than would be expected 
if the error bars were purely statistical. This is because 
the data at different times are obtained from the same 
experimental sequence, with the absorption integrated 
over different time windows, and so there is a high degree 
of correlation in the noise at different times. For the 
translational motion the best fit parameters are T c — 
21±1K, Tn = 329±117K, T i2 = 28±3K, n =84±14/is 
and T2 = 880 ± 160 /is. For the rotational motion, the fit 
parameters are less precise but are all consistent with the 
above values. 

We now explain why we choose to fit to this model, and 
we interpret the two time constants. In a simple ther- 
malization model [42 , the temperature of the molecules 
is given by 



dT 
~dt 



-R(T — Tb g )/ ft, 



(8) 



where R is the collision rate, Tb g is the buffer gas tem- 
perature and n = (M + m) 2 / (2Mm). The collision rate 
is given by the product of the helium density, nn e , the 
collision cross-section, <r, and the thermally-averaged rel- 
ative speed v, R = nn e crv. The helium atoms are much 
lighter than the YbF molecules and so, over the temper- 
ature range of interest here, v is dominated by the speed 
of the helium and has little dependence on the tempera- 
ture of the molecules. In this case, provided the collision 
cross-section depends only weakly on the temperature, 
the collision rate is approximately constant. Then the 



T = (T - T bg ) exp[-Rt/K] + T bg , 



(9) 



where To is the initial temperature of the molecules. This 
single exponential decay does not fit well to the data. 
Modifying this model to take into account the depen- 
dence of v on the molecule temperature does not improve 
the fit. 

This leads us to propose a more complicated model in 
which the molecules thermalise with the buffer gas as be- 
fore, but the buffer gas temperature is initially raised by 
the ablation process to Tb g ,o, an d then cools with a time 
constant r bg to the temperature of the cell walls, T c . The 
function that describes this model is the same as Eq. ^ , 
but with the buffer gas temperature now a function of 
time, T bg = (Tb g ,o — T c ). exp[— t/ r bg] + T c . This model 
fits the data very well. For the translational motion, this 
fit suggests that the cross-section for thermalization of 
the YbF with the helium is a = 7.1 ± 1.0 x 10" 20 m 2 , 
that the buffer gas is initially at Tb g ,o = 48±4K and 
that it cools to T c = 20±1 K with a time constant 
°f r bg = 915±180 /is. Fitting the same model to the 
rotational temperature data, we obtain similar fit pa- 
rameters: a = 6.8 ± 6.0 x 10- 20 m 2 , T bg , = 49±36K, 
T c = 19±1K, and r bg = 460±250/is. The energy re- 
quired to explain the 30K increase in buffer gas tem- 
perature corresponds to roughly 0.5% of the pulse en- 
ergy of the ablation laser being rapidly transferred to 
the helium buffer gas. A time constant of about 1 ms 
for re-thermalization of the buffer gas with the cold walls 
of the cell is as expected given the size of the cell and 
the thermal diffusivity of helium at these temperatures 
and pressures. However, an atom- molecule thermaliza- 
tion cross-section of about 7 x 10 _20 m 2 is smaller than 
we would expect in light of the measurements presented 



in Sec. IV It also seems unlikely that the time constant 
for rotational thermalization would be so similar to that 
of translational thermalization. Therefore, although the 
model fits to the data, it only does so with parameters 
that we find unrealistic. 

These observations suggest that the rotational and 
translational motions thermalize with the helium on a 
timescale that is shorter than the time resolution of the 
experiment, and the temperature evolution we observe in 
Fig. [lO] is entirely that of the helium. This implies that 
the heat of ablation deposited into the helium is dissi- 
pated with two different time constants, one being about 
10 times longer than the other. We suppose that the 
heat is initially deposited into a localized region close to 
the target. This hot helium then mixes with the cold 
helium so that, on the timescale Ti, the temperature be- 
comes roughly uniform throughout the cell, but higher 
than that of the cell walls. Finally, the entire bulk of 
helium cools to the cell temperature with a longer time 
constant r 2 . In terms of diffusion modes, the initial tem- 
perature distribution is represented by a sum of modes, 
the higher-order ones representing the more rapid varia- 
tions in temperature and decaying away rapidly, leaving a 
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FIG. 10: (Color online) Measured temperatures of the trans- 
lational motion (open points, red) and rotational motion 
(filled points, blue) when the buffer gas cell is maintained at 
approximately 20 K and the helium density is approximately 
0.9 x 10 22 m -3 . Some of the error bars are smaller than the 
size of the data points. The solid and dashed lines are double 
exponential fits to the translational and rotational data, as 
described in the text. 



small amplitude lowest-order mode which decays slowly. 
This model for the diffusion of heat is exactly analogous 
to the diffusion of the molecules in the cell, giving the 
double exponential decay that we have fit to the data, 
Eq. ([7]). Modelling the diffusion of heat to the cell walls 
using the thermal diffusivity of helium at this tempera- 
ture and density, we find that the time constants found 
from the fit are reasonable. We conclude that this is the 
most likely explanation for the cooling we observe. 



VI. MOLECULE DENSITY AND SATURATION 
OF THE ABSORPTION 

All of the measurements presented in this paper are 
based on absorption of a probe laser. In this section, and 
in Appendix [Aj we consider this absorption in more de- 
tail. We are particularly concerned to obtain an accurate 
value for the molecule density from the observed absorp- 
tion, and to understand how the fractional absorption 
depends on the probe intensity and on the collision rates 
between the molecules and the helium. 

^a) shows how the fractional absorption of the 
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probe laser depends on the laser intensity for three dif- 
ferent values of the buffer gas density. In these experi- 
ments, the temperature was 80 K. As the laser intensity 
increases, the fractional absorption decreases. This fall in 
fractional absorption is less rapid when the helium den- 
sity is high than when it is low. The data show that there 
is a saturation of the absorption, and that the intensity 
at which this saturation sets in increases with increasing 
buffer gas density. It is interesting to explore the mech- 
anism responsible for this saturation. 

In a two level system, the absorption saturates once 




0.00 



500 1000 1500 2000 2500 3000 3500 
Intensity at centre of probe laser (W/m 2 ) 




5 10 15 

Helium density (10 22 nr 3 ) 

FIG. 11: (Color online) (a) Fractional absorption versus probe 
intensity at three different helium densities. The temperature 
is 80 K. (b) Saturation intensity versus helium density at 80 K. 



the rate of excitation from ground to excited state ex- 
ceeds the rate of decay of the excited state. Once this 
occurs the rate at which photons are absorbed is limited 
by the decay rate, and an increase in the incident power 
does not produce a proportional increase in the absorbed 
power. Collisions increase the decay rate, resulting in 
an increase in the intensity required to reach saturation, 
accompanied by broadening of the spectral line. Our 
saturation data show a strong dependence on the helium 
density, but we observed no pressure broadening of the 
spectral lines in our saturated absorption spectra for den- 
sities up to 20 x 10 22 m -3 . We conclude that this is not 
the mechanism responsible for the saturation observed 
here. 

The molecules are not well modelled by a two level 
system, but they can be fairly well modelled using three 
levels. The probe laser connects levels 1 and 2, while 
level 3 represents all other levels in the molecule and is 
not coupled at all to the laser. Since level 2 can decay 
to level 3, the molecules are optically pumped into level 
3 by the probe laser. The intensity needed to saturate 
the absorption is set by a competition between the optical 
pumping of molecules out of level 1, and any process that 
puts population back into this level. We identify three 
such processes. Firstly, diffusion brings new molecules 
into the interaction region. Secondly, inelastic collisions 
transfer molecules from level 3 back to level 1. Thirdly, 
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when the Doppler width is much larger than the natural 
linewidth, optical pumping into level 3 occurs only in the 
velocity group that is resonant with the laser. Elastic 
collisions then repopulate that velocity group from the 
reservoir of molecules with all other velocities. The rate 
for the first process is inversely proportional to the helium 
density and is very small compared to all other relevant 
rates over the entire range of densities used here. The 
rates for the second and third processes are proportional 
to the helium density and, together, are responsible for 
the saturation of the absorption that we observe. 

The expected fractional absorption is calculated in Ap- 
pendix [A] To model our data, we first make the small ab- 
sorption approximation which simplifies the results con- 
siderably. In this approximation we assume that the 
laser-driven excitation rate is constant throughout the 
length of the absorbing column of gas, because the change 
in the intensity as the laser propagates through the col- 
umn is small. This allows us to replace the variable in- 
tensity / w ith th e constant value L m on the right hand 
side of Eq. flA24] ). The maximum fractional absorption 
is 30% for the data analyzed here and shown in Fig. 11 
This gives a rough indication of the maximum error that 
we make by assuming a constant excitation rate through- 
out the sample. Next we make the assumption, which 
we will verify later, that all the relevant collision rates 
are considerably smaller than the natural decay rate of 
the upper level, T. In this case, the fractional absorp- 



tion is given by Eq. ( A40 ) with the saturation intensity 
I s replaced by I' s which is given by Eq. (A41). The ratio 



of the Doppler width to the natural linewidth, w : enters 
the expression for the fractional absorption, and its value 
given by Eq. (A26) is w = 15 in this experiment. This be- 
ing much larger than 1, the fractional absorpti on va ries 
with intensity as (1 + ii n //g) -1 / 2 . Equation (A40) as- 



sumes a constant intensity across the laser beam but the 
probe beam used in these experiments had a Gaussian 
intensity distribution. We need to account for this in our 
analysis by integrating over the distribution. The result 
for the fractional absorbed power, which is the quantity 
measured in the experiment, is then 



A- 



f 

Jo 




pdp, (10) 



where Iq is the intensity at the centre of the probe beam. 
The quantity I r s is a modified saturation intensity and 
is given by Eq. (A41), while is the optical thickness 
of the sample for laser light resonant with the 1—^2 
transition and is given by Eq. ( A30 ) . We fit the data 



shown in Fig. 11 ^a), along with similar datasets taken 
at other helium densities, to Eq. (10) with ao and I' s as 
the only free parameters in the fit. The solid lines in 



Fig. 11 ^a) are the results of this fit for the three datasets 
shown there. 



Figure [TTJb) shows how the value of I' s obtained from 
these fits depends on the helium density, along with a 
weighted linear fit to these data. The linear dependence 
of the saturation intensity on the helium density supports 
our assumption that the collision rates are small com- 
pared with T, at least for the lower density data which has 
the greatest weight in the fit. When the collision rates 
become comparable with T, we expect terms of higher 
order in the density to appear in the expression for the 



saturation intensity, as shown by Eq. (A19). A weighted 
fit to quadratic and cubic polynomials shows the coeffi- 
cients of these higher-order terms to be zero within the 
error of the fit. The intercept of the linear fit is also zero 
within the fit error, as we would expect, and the values 
of ao are independent of helium density. We con clude 
that our model is the correct one. We see from Eq. ( A41 ) 



that the slope of the straight line fit gives us the ratio 
(713 + 7c)/^He, where 713 is the rate at which collisions 
transfer population from level 1 to level 3, and j c is the 
rate for velocity-changing collisions. We do not know 
which of these two rates, 713 or 7 C , is the larger, and our 
data does not help us to distinguish between them, so we 
simply take them together. Both rates are responsible 
for re-filling the hole produced by the laser in the com- 
bined velocity and internal-state distribution. We define 
the thermally-averaged cross-sections, cr 13 and cr c , via 



7x nax 



SkT 
M He 7r' 



(11) 



The sum of the two cross-sections is obtained from the 
gradient g = 1.0 x 10 -20 (Wm -2 )/m -3 of the straight 
line fit: 



&13 



3X 3 r(l-r)pg M Ue 7r 



2irhc 



8kT 



(12) 



Here, p is a factor that depends on the number of 
mag netic sub- states in levels 1 and 2, and is given by 
Eq. (A15). For this experiment its value is p = 1/3. The 



branching ratio r is the probability that the excited state 
decays back to the state resonant with the laser. For our 
experiment, it is the product of the Franck-Condon fac- 
tor, Z, for the — vibrational transition, and an angular 
factor that accounts for the probability of returning to 
the same rotational and hyperfine state. The calculation 
of this angular factor is similar to the one presented in 
the appendix of reference [43] . 

Consider driving a transition from a positive-parity 
ground state X 2 S + (v, AT, J = N ± 1/2) to a negative- 
parity excited state A 2 Hi/ 2 (v' , J' = N + 1/2), where 
the spin-orbit splitting of the excited state is large com- 
pared to the rotational splitting, there are no inter- 
mediate electronic levels, and we neglect for the mo- 
ment the hyperfine structure arising from any nuclear 
spin. The only possible decay routes are to the states 
X 2 £+(t/', N", J") where v" can take any value with 
weight given by the Franck-Condon factor for v' — » v n \ 
but (N" , J") are constrained to one of three possibilities, 
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(1) (N" = iV, J" = AT+1/2) (2) (TV" = iV, J" = N-l/2) 
or (3) (iV" = TV + 2, J" = TV + 3/2). The corresponding 
transitions are usually called Qn, R12, and P12, respec- 
tively. We find that the relative branching ratios between 
these three depends on the value of J'. When J' = 1/2 
the three branches have relative weights of 2/3,0,1/3, and 
as J' increases these tend to the asymptotic values 1/2, 
1/4, 1/4. For the data analyzed here the excitation was 
to J' = 27/2, and the ratios are 0.50, 0.24, 0.26, very 
close to their asymptotic values. In excitation, the two 
transitions with J = N ± 1/2, Qn and R42, are unre- 
solved in our Doppler-limited experiment, but are sepa- 
rated by more than the natural linewidth. Both transi- 
tions are excited, but the velocity group is different for 
the two excitation routes. In this case, we should take 
the mean value of r for the two branches. The situation 
is further complicated by the fact that each transition 
consists of a pair of hyperfine components which are un- 
resolved in the Doppler-limited spectrum. In the case of 
the Qn line the components are resolved in the Doppler- 
free spectrum and we should take the two hyperfine states 
as distinct. This halves the hyperfine- averaged value of 
r for this line. For the R12 line, the situation is different 
because the hyperfine components are unresolved in the 
Doppler-free spectrum. In this case it is better to treat 
the two hyperfine levels as a single level and the value of 
r is unchanged. These considerations lead us to an aver- 
aged branching ratio of r « 0.25Z. We do not know of 
any direct measurements or calculations of the Franck- 
Condon factor for the A(V = 0)-X(V' = 0) transition of 
YbF. However, we know that when the molecule is ex- 
cited strongly on the F = 1 component of the Qn(l/2) 
line, the mean number of photons scattered is 1.9 ± 0.1 
[44] . This information, together with the table of branch- 
ing ratios given in [43] provides an estimated value of 
Z = 0.91 ±0.05. A high value for Z is to be expected 
since YbF closely resembles the alkaline-earth monofluo- 
rides which all have this property [45 . Finally then, we 
take r ~ 0.23 for this experiment. 



Putting these values into Eq. ( 12 ) we obtain a is -fa, 



37 x 10- 20 m 2 . There are a number of uncertainties in 
determining this value. The fractional error on the gra- 
dient of the straight line fit is 15%. Using the spread 
in the diffusion lifetimes measured at the same nominal 
helium density we estimate the fractional uncertainty in 
our measurement of the helium density to be 20%. The 
fractional uncertainty in calibrating the laser intensity is 
also 20%. The fractional uncertainty in the branching 
ratio r is estimated to be 24%, partly from the 5% un- 
certainty in the Franck-Condon factor, but mostly from 
the uncertainty of how best to treat the four unresolved 
components of the line as discussed above. This leads to 
a 19% error in the cross-section which depends on the 
value of r(l — r). Other possible errors come from the 
approximations we have made in our model. First, we 
have assumed that the absorption is small; this does not 
contribute a significant error at this level. Second, our 
analysis relies on the approximation that the collision 



rates are much smaller than T. We have previously esti- 
mated T to be approximately 6 x 10 T rads _1 [22]. When 
we use our cross-section to calculate the collision rate, 
we obtain 713 + 7 C = 5 x 10 T s _1 at n = 2 x 10 23 m -3 . 
This implies that the collision rate is approaching T at 
the upper end of the plot. However, it is the points at 
low helium densities that primarily fix the gradient of 
the fit, since these have far more weight in the fit, and 
for these values the collision rates are much smaller than 
T, consistent with our original assumption. Third, we 
have used a fitting function that is appropriate when the 
laser is resonant for molecules at rest, but this cannot be 
satisfied for all four of the unresolved components. We 
find that setting 5' L = w/2 decreases the cross-section by 
3%, and setting 5' L = w decreases the cross- section by 
12%. Considering the mean displacement of the compo- 
nents from the line centre at 80 K, we conclude that the 
resulting change in the cross-section, and the associated 
error, are both negligible for this measurement. The final 
result is <ti 3 + a c = 37 ± 14 x 10" 20 m 2 . This result is 
close to our calculated YbF-He diffusion cross-section at 
80 K, given in Table [I] as one would expect if velocity- 
changing collisions dominate over inelastic collisions in 
determining the saturation intensity. 

The values of ao obtained from the fits to the data 
shown in Fig.JTTJa), along with similar datasets taken 
at other helium densities, are independent of helium 
density and are consistent with one another across the 
whole range of helium densities. The mean value is 
ao = 18.1 ± 0.8. The fact that for the four unre- 

solved components has a much larger effect on ao than 
on l£, and by considering the effect of this we obtain a 
corrected value of ao = finer L = 23 ±4. Using the above 
values of r and p, the val ue fo r the optical absorption 
cross-section given by Eq. (A14) is a = 1.1 x 10 _14 m 2 . 
Taking an absorption length of L = 3 cm and the frac- 
tion of all molecules in iV = 13 to be fi = 0.05 at 
80 K, we find the total density of YbF molecules to be 
n ~ 1.3 x 10 18 m -3 in these experiments. The total num- 
ber of YbF molecules produced is approximately 3 x 10 13 . 
Of course, the YbF density depends on the energy and 
spot size of the ablation laser, the target quality and the 
buffer gas density. Nevertheless, we consistently produce 
this density of molecules once the experimental condi- 
tions are optimized. 



VII. CONCLUSIONS 

We have investigated the dynamics of YbF molecules 
produced inside a helium buffer gas cell using a com- 
bination of absorption imaging, absorption spectroscopy 
and saturated absorption spectroscopy. While we have 
focussed on this particular molecule, we expect our re- 
sults to carry over to a wide range of other polar radicals 
produced using similar methods. 

The absorption images show that, upon ablation of the 
target, the molecules that are formed expand ballistically 
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into the cell on a /is timescale. When the helium den- 
sity is low, this initial ballistic expansion completely fills 
the cell, but at higher helium density the expansion is 
arrested with the molecules confined to a region around 
the target. The molecules then diffuse slowly to the walls 
of the cell. We have measured how the diffusion time de- 
pends on the helium density and find the expected linear 
dependence at low density. At higher densities the dif- 
fusion times no longer increase with density but instead 
reach a roughly constant value. Our modelling of the ex- 
periment shows that the diffusion is more rapid when the 
molecules are localized near the target than when they 
fill the cell, and that in the former case the decay curves 
depend strongly on the probe position, with the density 
decaying much faster when the probe is close to the tar- 
get. This observation, along with the insights obtained 
from the absorption imaging, suggest that the levelling 
off of the lifetimes is due, at least in part, to the initial 
confinement of the molecules to a region close to the tar- 
get. By comparing the diffusion of YbF and Li at low 
helium densities, and using values for the Li diffusion 
cross-section calculated from a well-established interac- 
tion potential, we have measured the YbF-He diffusion 
cross-section at two temperatures. The experimental val- 
ues agree with the theoretical values we calculate from 
the YbF-He potential given in [40]. This is a first ex- 
perimental test of that potential. We have compared our 
experimental results for the Li diffusion times with the 
results of a numerical model of diffusion in this cell. The 
results agree at 80 K and 20 K, but differ by a factor 
of two at 293 K where the results are most precise. This 
discrepancy is not understood but suggests that diffusion 
cross-sections determined by measuring lifetimes in a cell 
may not be accurate unless a reference atom is used. 

We have measured how the translational and rotational 
temperatures of the molecules change as a function of the 
time since their formation, for typical helium densities of 
about 10 22 m -3 , and have compared our experimental 
results with those of various thermalization models. Our 
conclusion is that the YbF molecules thermalize with the 
helium more rapidly than the 50 /is time resolution of 
the experiment. Instead of observing that thermalization 
process, we see the cooling of the helium buffer gas which 
heats up when the target is ablated and then slowly cools 
back down to the temperature of the cell walls. We see 
that the cooling occurs on two separate timescales which 
differ by a factor of 10, and we suppose that the heat of 
ablation is deposited in a local region close to the target, 
and that this heat first dissipates quickly into the rest of 
the gas, and then more slowly as it diffuses to the walls. 

Finally, we have made a detailed study of the absorp- 
tion of a probe laser by molecules in a buffer gas cell, 
elucidating the main dynamical processes that are im- 
portant. As the power of the probe laser is increased 
the absorbed power increases at first, but eventually sat- 
urates because the population in the resonant state is 
optically pumped into a multitude of other states that 
are not resonant with the laser. The resulting hole in the 



combined velocity and internal state distribution is re- 
populated by collisions, and so saturation of the absorp- 
tion sets in when the rate of optical pumping is equal to 
the relevant collision rates. We have measured the frac- 
tional absorption as a function of laser intensity over a 
range of helium densities and find that the data fits well 
to our model. The intensity needed to saturate the ab- 
sorption is directly proportional to the helium density for 
densities up to 10 22 m -3 . The cross-section for the "hole- 
filling" collisions is consistent with our theoretical value 
for the YbF-He diffusion cross-section, as we would ex- 
pect if velocity- change is the dominant hole filling mech- 
anism. The same data are used to estimate the total 
YbF density, which is approximately 10 18 m -3 . The to- 
tal number of YbF molecules produced is approximately 
3 x 10 13 . 

Using an improved setup, we have now reduced the 
cell temperature to 3.5 K. At this temperature there will 
be about 3 x 10 12 molecules in the rotational ground 
state. Using a hydrodynamic flow of helium to sweep 
the molecules out of the cell, we expect that at least 
10% can be extracted into a beam with a divergence of 
about O.lsteradian [3T] . resulting in a cold, slow- moving 
molecular beam with a flux of 3 x 10 12 ground state 
YbF molecules per steradian per pulse. This is 2000 
times more intense than the flux produced in a super- 
sonic source of the same molecules [44], and provides 
molecules at considerably lower speed, and so promises a 
large improvement in the sensitivity of molecular beam 
experiments such as the measurement of the electron's 
electric dipole moment. 



Acknowledgments 

We are grateful to Jon Dyne, Steve Maine and Va- 
lerijus Gerulis for their expert technical assistance. We 
thank Roman Krems for sending us the YbF-He inter- 
action potential. This work was supported by the EP- 
SRC, the STFC, and the Royal Society. The research 
leading to these results has received funding from the 
European Community's Seventh Framework Programme 
FP7/2007-2013 under the grant agreement 216774. 



Appendix A: Laser absorption by molecules in the 
presence of a buffer gas 

To understand the interaction of the probe laser with 
the molecules in the buffer gas, we use a 3 level model of 
the molecule. In this model, a plane polarized laser field, 
E = Eoz cos(cjt), is detuned from resonance with the 
transition from level 1 to level 2 by an angular frequency 
5. Level 3 represents all other relevant molecular levels 
and is not coupled to the laser field. Level 2 can decay 
by spontaneous emission, either to level 1 at rate rT, 
or to level 3 at rate (1 — r)T. Inelastic collisions with 
helium atoms also cause level 2 molecules to decay to 
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levels 1 and 3, with rates 721 and 723 respectively. Elastic 
collisions dephase level 2 at a rate j e i. Levels 1 and 3 
are assumed to be in the ground electronic state of the 
molecule, so spontaneous emission out of these levels is 
negligible. Collisions result in population transfer from 
level 1 to level 3, and vice versa, with rates 713 and 731 
respectively. There is not enough energy for a collision to 
excite molecules to level 2. Molecules move in and out of 
the interaction region with a characteristic rate 0, which 
is the inverse of the average time taken to diffuse out of 
the probe laser beam. Following the usual treatment that 
leads to the optical Bloch equations for a two level atom, 
modified to account for the third level, for the collisions, 
and for the flux, we obtain the following equations for 
the rate of change of the level 1 and 2 populations, p\\ 
and P22, and of the coherence between these two levels 

P\2- 



P11 = ^ (P12 ~ P21) + (rT + 721) P22 ~ 7i3Pii 

+73i (1 - P11 - P22) - #11 + A0, (Al) 

P22 = -i^ (P12 - P21) - (r + 721 + 723) P22 - #22, 



Pl2=i^ (Pll-P22)+M*-^ (r + 721+723)-7el ) Pl2- 



(A2) 

'12- 

(A3) 



Here, f\ is the fraction of all molecules in level 1 when 
the probe laser is off, and Q is the Rabi frequency defined 
by Q = E0Z12/T1 where z±2 is the z component of d\i — 
(l|d|2) and d is the dipole moment operator. As the laser 
propagates a small distance Az through the molecular 
sample of number density n, its intensity, / = ccqEq/2, 
changes due to absorption and stimulated emission by 
the amount 



AI = ihuj^- (pi2 - P21) nAz. (A4) 

It will be useful to define the dimensionless rates 
i i} = jij/T, i el = lel /T, 6' = S/r, cj>' = cj>/T and I*, = 
1 + 72i + 723 + ^7e^' a l° n g with a dimensionless quantity 
proportional to the laser intensity, s = 2Q 2 /T 2 = I/I S1 
where I s is called the saturation intensity. It will also 
be helpful to make use of the relationship between 731 
and 713 which is easily found by ensuring that the pop- 
ulations are in a steady state when the laser is turned 
off: 



713 
731 



1-/1 

h 



(A5) 



We will calculate the absorption averaged over a time 
window that is very short compared to the characteristic 
diffusion time, but very long compared to 1/r. In this 
case, we can take n to be a constant, and we can use the 
steady state solution to the above dynamical equations. 



After some algebra we find the steady state values for the 
density of molecules in levels 1 and 2 to be 



n 2,. 



(r + 721 + 723 -+ 

R 



(A6) 



x = (713 + 731 + <t>) (r + 721 + 723 + </>) , 
Y = R (713 + 723 + 2731 + T(l - r) + 20) , (A9) 



(A7) 
(A8) 



where R is the laser excitation rate, 



R 



V 

1 rri 



2 V 4<5' 2 + 



(A10) 



With the populations in steady-state, the equation for 
the intensity reduces to 



dl 

dz 



-huj(T + 721 + 723 + 0)n 2 , a 



(All) 



It is convenient to rewrite this equation in terms of the 
reduced intensity s = 2^ 2 /r 2 , which we can do using the 
relationships between /, Eo, and ft given above. We also 
introduce a reduced length Z — z/L, L being the total 
length of the interaction region. Thus we obtain 



ds 
dZ 



-2aL(l + 7 21 +7 23 + <//)n 2 , ss , 



where 



3A 2 



(7 = 



(A12) 



(A13) 



is the absorption cross-section of the molecule. To reach 
Eq. (A13) we have used the expression for the decay rate, 
r = siveohc 3 ^k^k^jk^ have made the approximation 
that the frequencies of all allowed transitions are equal, 
ujjk = 27rc/A, and have generalized the notation in order 
to handle the prickly issue of degeneracy. To represent 
the molecule accurately, we have to allow each of levels 1 
and 2 to repres ent a set of degenerate or near-degenerate 
levels. In Eq. (A13), i represents a particular Zeeman 
sub- level of level 1, i = (1, M p ), j is the sub-level (2, M p ) 
of level 2, and k is an index which runs over all the sub- 
levels of levels 1 and 3. The branching ratio r is the 
rate for spontaneous decay to any of the sub-levels in 1, 
divided by the t otal s pontaneous decay rate. Using this, 
we see that Eq. (A13) can be written in the form 



where 



P 



3X 2 rp 



\(l,M p \d.z\2,M p )\< 
Emi l(2jMp|d|l,Mi] 



(A14) 



(A15) 



is a factor that accounts for the degeneracy, the sum run- 
ning over all sublevels Mi of level 1. Finally, when level 
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1 is degenerate and the molecular sample is unpolarized, 
we have to remember to take the average value of a cal- 
culated for each of the ground state sub-levels M p . Using 
these same definitions, we can write the saturation inten- 
sity as 



nhcT 
3X 3 rp ' 



(A16) 



R etur ning to the calculation of the absorp tion we use 
Eq. (A7) to substitute for n2, ss in Eq. (A12), leading us 



to the result 



ds as 
dZ = ~T+05' 



where 



a = 



final 



r 2 



p/2 
1 T 



r 2 + 45' 2 



(A18) 



and 



a = IT r + ^3 + 723 + 2731 + 20' 



2r^( 7 ; 3 + 7 3 i + 0O(i + 72i + 723 + 00 V^ 2 + 4£' 2 



1 T 



00 



r /2 

1 T 



+ 4<5' 2 



(A19) 



These equations define a^ and /3q which are the values of 
a and f3 when J = 0. 



Integrating Eq. (A17) we obtain the result 



In [ — ] + /3(s - s in ) 



(A20) 



where S[ n is the value of the reduced intensity at the 
entrance of the cell. We have taken n (and hence a) to be 
independent of Z over the interaction length. If it is not, 
we have simply to replace the product nZ by the column 
density J n(Z)dZ. We can rewrite this transcendental 
equation in the form 



Pse p * = /3s 



-aZ+/3s in 



which has the solution 



s(Z) = ^w(j3s- m e 



(-aZ+(3s in ) 



(A21) 



(A22) 



Here W(x) is the Lambert W function 2 . The fractional 
absorption measured at the end of the interaction region 
(i.e. at Z = 1 ) is 



^in -^out ^ 

S in PS 



— W((3s in e 

Sin V 



(-a+0Sin) 



). (A23) 



2 The Lambert W function is sometimes known as the omega func- 
tion or the product log. It is the inverse function of f(W) = 
We w 



So far we have neglected the distribution of speeds in 
the gas. Therefore, our results are only valid in the case 
where the Doppler width is small compared to T. This 
is not often the case, and so we need to extend the treat- 
ment to include the distribution of speeds, g(v), where v 
is the component of velocity along the propagation direc- 
tion of the probe laser. This modifies the treatment in 
two ways. First, we have to write the equations for the 
absorption as a function of v and then integrate over v at 
the end to obtain the total absorption. Let us use a new 
symbol, 5l, for the detuning of the laser from the at-rest 
transition frequency, and use 5d for the Doppler shift of 
a molecule, = 2ttv/X. Then, S = Sl — Sp. Instead of 



(A17) Eq. flA17D, we then have 



ds 
dZ 



a(v)g(v)s 
1 + P(v)s 



dv, 



(A24) 



where a( v) an d /3(v) are the same as given by Eq. flA18l ) 
and Eq. ( A19| ) but with 5' replaced by 5' L 
other symbols have the same meaning as before, 
velocity distribution is 



5' D . All the 
The 



g{v) dv 



M 



2irkT 
1 



exp 



-Mv 2 
2kT 



dv 



exp 1 -t§ ) d6 f D . 



(A25) 



Here, the width w is related to the temperature T and 
the molecular mass M through 



w 



8ir 2 kT 
M\ 2 Y 2 ' 



(A26) 



The second modification that we need to make is to 
account for velocity-changing collisions in the rate equa- 
tions. The probe laser tends to modify the velocity distri- 
bution by depleting the resonant velocity group, whereas 
collisions tend to restore the velocity distribution towards 
thermal equilibrium. The extra terms that appear in the 
rate equation for n\ are 

l n i( v )] c = - n i( v ) J lcW,v)dv f + J 'y c (v,v')ni(v')dv', 

(A27) 

where 7 C (^ / , v)dv' is the rate for collisions that change 
the velocity from v to a value in the range v r — >• v r + dv' . 
We make the "strong collisions" approximation, meaning 
that 7 C (v' , v) is independent of v [46]. Then, to make sure 
that the system is in equilibrium when the probe laser is 
turned off, we must have 

- ni (v) J j c (v')dv' + lc {v) J ni [W)dW = 0, (A28) 

from which it follows that j c ( v ) — g( v )la where 7 C is 
the total collision rate integrated over all final velocities, 
7 C = J j c (v f )dv f . It follows that the extra terms that 
appear in the rate equation for n\ reduce to 



l n i( v )] c = ~ n i( v hc + fig(v)j c n. 



(A29) 
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If we look back at Eq. (Al), re- write it as a function 



of v and multiply through by n so that it becomes a rate 
equation for ni(v), we find that the last two terms in the 
equation are —n\(v)(j) + fig(v)(j)n, which has exactly the 
same form as the above equation. A similar replacement 
occurs in the equation for ^22- So we see that to include 
velocity-changing collisions we simply have to replace (j) 
with (j) + 7c. 

To summarize, the change in the probe intensity is 
given by Eq. (A24) where g(v) is given by Eq. (A25), 



a(v) = a T'S/lT^+4S'(v) 2 ), p(v) = A,r£/(r£+4^) 
and ao and /3q are 



a 



finer L 

L rp 



(A30) 



l-r + 7i 3 + 723 + 273i+27c + 20' 



2r^( 7 i 3 + 731 + l'c + 00(1 + 721 + 7 23 + l'c + </>0 ' 



These are very general results that can be used to cal- 
culate the degree of absorption for almost any case, once 
the relevant rates are known. The same equations give 
the spectral lineshapes to be expected as the laser fre- 
quency is scanned, and tell us how the fractional ab- 
sorption changes with the intensity of the probe laser. 
In general, we will have to solve Eq. ( |A24 ) numerically. 
Furthermore, if the intensity of the probe laser is not uni- 
form over the interaction volume, we should divide up the 
volume into small elements and integrate over these too. 
However, in many cases of interest the results simplify 
considerably because some particular quantity is small. 
Next, we examine some of these simpler cases. 



1. The low intensity limit 

Whe n Psj n < 1 we can neglect it in the denominator 
of Eq. ( A24 ) , so that we have 



= — f a(v)g(v)sdv = —as 
dZ J 

where 

a = J a(v)g(v) dv. 
The fractional absorption is 

'-'out/ Sh 



(A32) 



(A33) 



A=l 



1 



The full expression for a is quite complicated, but when 
the laser is on resonance with the atomic transition (i.e. 
S = Sd) it simplifies to a reasonably compact form: 



a 



a 



^e r ?/K)[l -erf(rV(2w))]. (A35) 



2w 



When w ^> T' T , which will usually be the case if the 
gas is warm and the collision rates are not too much 
larger than T, this expression simplifies further to a = 
a ^/7rT / T /(2w). 



2. The low absorption limit 

When the absorption is small the value of s hardly 
cha nges a nd it is valid to set s = S[ n in the integral of 
Eq. (A24). Then, the fractional absorption is 



A 1 "^out/^in 



1 + P(v)s iT 



dv. 



(A36) 



When the laser is on resonance, this integral evaluates to 
an analytical form, 



A 



r(l + /3 Sin) 



2w VI + /?0«ii 



1 rp 



1-erf [ ^Vl + A)Si, 
2w 



(A37) 

When the temperature is high so that w 3> T' T , and 
the intensity is not too high, the expression simplifies to 



(A31) A = ao ^lV(2uVl + A)Sin). 



3. The two level limit with no collisions and no flux 

To verify that we can recover some familiar results, let 
us specialize to the case of just two levels by setting all 
the rates involving level 3 to zero and by setting r = 1. 
Neglecting all the remaining collision rates and the flux 
term 0, we obtain 



finer L 
1 + 4£' 2 ' 
1 

1 + 4(5 /2 ' 



When the intensity is low and the velocity distribution is 
neglected, we have, using Eq. (A34) 



A=l — exp 



finer L\ 
+ 4£' 2 J 



(A38) 



When the absorption is small and the velocity distribu- 
tion is neglected, but the input intensity is not necessarily 
small, we have 



A = 



finer L 



l + 4J' 2 +/ in // s 



(A39) 



where I s is the saturation intensity defined by Eq. (A16). 



( A34) I n Eqs. ( A38 ) and ( A39 ) we have kept f\ in the numerator 



even though it is unity in this case because, as we shall 
soon see, the same equations hold for other cases. We 
note that in this two level limit, the absorption cross- 
section may reach its maximum possible value, a max = 
3A 2 / (27r). For example, in an atomic transition with zero 
angular momentum in the ground state and one unit in 
the excited state (an s-p transition) we have p = 1 and 
so er = a max . By contrast a p-s transition has p = 1/3 
for the M = component and p = for the other two 
components meaning that for an unpolarized sample the 
level- aver aged absorption cross-section is er = <r max /9. 
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When we include the velocity distribution and put the 
laser on resonance, we have 



A 



a/tt hnaL 



1+ 



2w 



Jin 
I, 




(A40) 



4. The low pressure, very low flux limit 

Returning to the 3 level system, let us assume that all 
the collision rates are much smaller than (1 — r)T, and let 
us also take the flux to be much smaller than the colli- 
sion rates. This limit applies to most of the experiments 
reported in this paper. 

In almost all situations we expect because the 

effective degeneracy of state 3 is vastly larger than that 
of level 1. Then we have 713 ^> 731. The only exception 
is when level 1 is the ground rotational state and the 
rotational temperature is much lower than the rotational 
spacing. Using this approximation, along with the ones 
discussed above, a and j3 simplify to 



a 



A) 



1 — r 

2(7i 3 +y c )" 



When the laser intensity is low and we neglect the dis- 
tribution of velocities the fractional absorption is again 
given by Eq. (A38), just as in the two- level case. Simi- 
larly, when the absorption is small and we neglect the ve- 
locity distribution, we again have Eq. (A39) except that 
I s in this equation is replaced by a new saturation inten- 
sity 



1: 



2 713 + 7c 



1 



I, 



(A41) 



When the absorption is small and we include t he ve locity 
distribution, the result is the same as i n Eq . (A40), but 
with I s again replaced by the I' s of Eq. ( A41| ). 



We see that for the three level system in this limit, 
the intensity required to saturate the absorption is pro- 
portional to the sum of two rates, the collision rate that 
transfers population between levels 1 and 3, a nd the rate 
for ve locity changing collisions. Equations (A10) and 
(A41) together show that, when 6 = and / = I' s1 we 
have (1 — r)R = 713 +7 C (in this low pressure limit where 



= 1). Since 1/(1 — r) is the mean number of photons 
scattered by the molecule before being pumped into level 
3, (1 — r)R is the optical pumping rate. So I' s is the in- 
tensity required to make the optical pumping rate equal 
to the rate at which collisions redistribute population be- 
tween levels 1 and 3 or between one velocity group and 
another. 

5. The low flux, very low pressure limit 

In this limit we take <j) <C (1 — r)Y and take all the 
collision rates to be far smaller than <p. In this case we 
have 

a « finer L, 
1-r 

Po 



Icb' 

Once again, the fractional absorption in both the low 
intensity limit and the low absorption limit is the same 
as in the two level case, but with I s now replaced by 



I" 



rT 



J, 



(A42) 



The intensity required to saturate the absorption is pro- 
portional to (j). When 6 = and / = I'J, we have 
(1 — r)R = 0, showing that I'J is the intensity for which 
the optical pumping rate equals the rate at which new 
molecules enter the interaction region. 



6. Summary 

When the intensity is low, the limiting rate is the 
laser excitation rate and the fractional absorption is 
A = 1 — e~ a . The result does not depend on flux, 
velocity-changing collisions, rotation-changing collisions 
or vibration-changing collisions. Collisions only enter 
through the modification of the decay rate of level 2 from 
r to IV. When the intensity is not low the result is more 
complicated, but it simplifies if the degree of absorption 
is low. In that case the fractional absorption on reso- 
nance has the form given by Eq. ( A40 ) where only the 



expression for the intensity needed to saturate the ab- 
sorption changes from one limiting case to another. This 
saturation intensity is proportional to T in the two-level 
limit, to 713 + 7 C in the low pressure, very low flux limit, 
and to cb in the low flux, very low pressure limit. 
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